Elastic properties of hydrogenated graphene 
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There exist three conformers of hydrogenated graphene, referred to as chair-, boat-, or washboard- 
graphane. These systems have a perfect two-dimensional periodicity mapped onto the graphene 
scaffold, but they are characterized by a sp 3 orbital hybridization, have different crystal symmetry, 
and otherwise behave upon loading. By first principles calculations we determine their structural 
and phonon properties, as well as we establish their relative stability. Through continuum elasticity 
we define a simulation protocol addressed to measure by a computer experiment their linear and non- 
linear elastic moduli and we actually compute them by first principles. We argue that all graphane 
conformers respond to any arbitrarily-oriented extention with a much smaller lateral contraction 
than the one calculated for graphene. Furthermore, we provide evidence that boat-graphane has a 
small and negative Poisson ratio along the armchair and zigzag principal directions of the carbon 
honeycomb lattice (axially auxetic elastic behavior) . Moreover, we show that chair-graphane admits 
both softening and hardening hyperelasticity, depending on the direction of applied load. 

PACS numbers: 81.05.ue, 62.25.-g, 71.15.Nc 
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I. INTRODUCTION 

The hydrogenated form of graphene is referred to as 
graphane. It is described as a two-dimensional, periodic, 
and covalently bonded hydrocarbon with a C:H ratio of 1. 
Hydrogen atoms decorate the carbon honeycomb lattice 
on both the top and bottom side (see Fig.[T]). Graphane 
was theoretically predicted by Sofo et a/., 1 further in- 
vestigated by Boukhvalov et al£ and eventually grown 
by Elias et al£ The investigation of graphane properties 
was originally motivated by the search for novel materials 
with possibly large impact in nanotechnology. 

The attractive feature of graphane is that by variously 
decorating the graphene atomic scaffold with hydrogen 
atoms (still preserving periodicity) it is in fact possi- 
ble to generate a set of two dimensional materials with 
new physico-chemical properties. This is obviously due 
to change in the orbital hybridization which, because of 
hydrogcnation, is now sp 3 -like. For instance, it has been 
calculated^ that graphane is an insulator, with an en- 
ergy gap as large as ~ 3 eV, while graphene is a highly 
conductive semi-metal. In case the hydrogenated sample 
is disordered, the resulting electronic and phonon proper- 
ties are yet again different 4 Hydrogenation likely affects 
the elastic properties as well. Topsakal et alA indeed 
calculated that the in-plane stiffness and Poisson ratio 
of graphane are smaller than those of graphene. In ad- 
dition, the value of the yield strain is predicted to vary 
upon temperature and stoichiometry. 

As far as the mechanical properties of graphane are 
concerned, the sp 2 -to-sp 3 change in orbital hybridiza- 
tion causes a major difference with respect to graphene. 
There in fact exist graphane conformers which are not 
isotropic, at variance with graphene which is so (in lin- 
ear approximation 5 ). This feature stimulates an intrigu- 



ing change of perspective, namely: hydrogenation could 
not only affect the actual value of some linear elastic 
moduli;* it could even dramatically change the overall 
mechanical behavior of the system by introducing an 
anisotropic dependence of its respose to an external load. 
This is in fact what we predict in this work by first prin- 
ciples total energy calculations, combined to continuum 
elasticity: we show that there is a graphane conformer 




Figure 1: (color online) Pictorial representations of the 
graphane conformers, obtained by different hydrogen deco- 
rations (the actual atomic positions are reported in Fig. [2]| . 
Top hydrogen atoms are indicated by red (dark) circles, while 
bottom ones by gray (light) circles. Shaded areas represent 
the unit cell and the corresponding lattice vectors are indi- 
cated by a and b. Panel a: graphene scaffold (full lines) with 
zigzag (zz) and armchair (ac) directions. Panel b, c, and d: 
chair-, boat-, and washboard-graphane, respectively. 
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(i.e., boat graphane as detailed below) showing a vanish- 
ingly small (possibly negative) Poisson ratio upon load- 
ing along given directions. In other words, we provide 
evidence that upon suitable hydrogenation a graphene 
sheet behaves as an axially auxetic material^ namely: it 
does not shrink, but actually slightly elongates perpen- 
dicularly to an applied traction force. Nonlinear elastic 
features show an interesting anisotropic behavior as well. 

This paper is organized as follows. In Sec. [Til the meth- 
ods and the general computational setup adopted in our 
calculations are outlined. In Sec. IIIII we provide a full 
structural characterization of three graphane conformers 
and we discuss their stability. In Sec. |TV] and Sec. [V] 
we describe their linear and nonlinear elastic properties, 
respectively, and we compute all the relevant elastic mod- 
uli. 



(a) C- graphane 




(b)B-graphane 

hi z 




(c)W-graphane 




Figure 2: (color online) Perspective representations of fully re- 
laxed graphane conformers. Gray (light gray) and red (dark 
gray) spheres represent carbon and hydrogen atoms, respec- 
tively. Labels C„ and H n (with n — 1,2,3 and 4) provide 
the atom identifications used in Table [I] Right panels show 
the orientation with respect to the armchair (ac) and zigzag 
(zz) direction, as well as the structural parameters h and ip 
reported in Table U 



Table I: Space groups and structural parameters for each 
graphane conformers. The cell parameters a and b are de- 
fined in Fig. Q] while the other quantities are reported in Fig.H 
Note that the B-graphane shows two types of C-C bonds while 
W-graphane exhibits a large buckling parameter, h. 
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II. COMPUTATIONAL SETUP 

All calculations have been performed by Density Func- 
tional Theory (DFT) as implemented in the Quantum 
ESPRESSO packaged The exchange correlation poten- 
tial was evaluated through the generalized gradient ap- 
proximation (GGA), using the Vanderbilt ultrasoft pseu- 
dopotential PW91. 8 A plane wave basis set with kinetic 
energy cutoff as high as 50 Ry was used and in most cal- 
culations the Brillouin zone (BZ) has been sampled by 
means of a (18x18x3) Monkhorst-Pack grid. The atomic 
positions of the investigated samples have been optimized 
by using the quasi-Newton algorithm and periodically- 
repeated simulation cells. Accordingly, the interactions 
between adjacent atomic sheets in the supercell geometry 
was hindered by a large spacing greater than 10 A. 

The elastic moduli of the structures under considera- 
tion have been obtained from the energy- vs-strain curves, 
corresponding to suitable sets of deformations applied to 
a single unit cell sample. As discussed in more detail 
in Sec. IIVI and Sec. El for any deformation the magni- 
tude of the strain is represented by a single parameter (. 
The curves have been carefully generated by increasing 
the magnitude of C in steps of 0.001 up to a maximum 
strain ICmairl = 0.05. All results have been confirmed by 
checking the stability of the estimated elastic moduli over 
several fitting ranges. The reliability of the above com- 
putational set up is proved by the estimated values for 
the Young modulus and the Poisson ratio of graphene, 
respectively 344 Nm" 1 and 0.169, which are in excellent 
agreement with recent literature^— Similarly, our re- 
sults for the same moduli in C-graphane (respectively, 
246 Nm' 1 and 0.08) agree very well with data reported 
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in Ref. H 

The stability of the three graphane conformers 
has been established by calculating the corresponding 
phonon dispersions. Phonon dispersions, have been ob- 
tained by means of Density- Functional Perturbation The- 
ory (DFPT);i2 based on the (2n + 1) theorem. In this 
case, during the self-consistent field calculation, the BZ 
has been sampled by a (16x16x3) Monkhorst-Pack grid. 
The accuracy of the phonon dispersion evaluations has 
been tested on a graphene sample (see below). 



III. STRUCTURE AND STABILITY OF 
GRAPHANE CONFORMERS 

By hydrogenating a honeycomb graphene lattice, three 
ordered graphane structures can be generated, namely: 
the chair (C-graphane), boat (B-graphane) and wash- 
board (W-graphane) conformer o 1 ' 13 shown in Fig. [5] 

Each conformer is characterized by a specific hydro- 
gen sublattice and by a different buckling of the carbon 
sublattice. In particular: in C-graphane the hydrogen 
atoms alternate on both sides of the carbon sheet; in B- 
graphane pairs of H- atoms alternate along the armchair 
direction of the carbon sheet; finally, in W-graphane dou- 
ble rows of hydrogen atoms, aligned along the zigzag di- 
rection of the carbon sublattice, alternate on both sides 
of the carbon sheet. A perspective view of the conformers 
is shown in Fig. [5] and the corresponding structural data 
are given in Table HI In C-graphane and W-graphane the 
calculated C-C bond length of 1.54 A is similar to the 
sp 3 bond length in diamond and much larger than in 
graphene. Moreover, we note that the B-graphane shows 
two types of C-C bonds, namely: those connecting two 
carbon atoms bonded to hydrogen atoms either lying on 
opposite sides (bond length 1.57 A) or lying on the same 
side of honeycomb scaffold (bond length 1.54 A). Finally, 
the C-H bond length of 1.1 A is similar in all conformers 
and it is typical of any hydrocarbon. 

The stability of the three graphane conformers has 
been established by calculating the phonon dispersion 
curves reported in Fig. [3] Graphene phonon spectrum 
is reported as well for comparison. No soft modes (with 
negative frequency) corresponding to possible instabili- 
ties were found along any high-symmetry direction of the 
Brillouin zone. Furthermore, as expected^ the zone- 
center longitudinal (LA) and transverse (TA) acoustic 
branches show a linear dependence upon the wavevector, 
while the acoustic mode ZA (with displacement patterns 
along the z-direction shown in Fig. [5]) shows a quadratic 
dependence. We observe that in C-graphane, as well as in 
graphene, the speed of sound (i.e. the slope of the acous- 
tic branches at T-point) is the same along the T — M 
and r — K directions. On the other hand, the B- and W- 
graphane conformers are characterized by different sound 
velocities along the T—X and T—Y directions. This is the 
fingerprint of an unlike elastic behavior: as extensively 
discussed in Sec. IIV1 C-graphane is elastically isotropic 
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Figure 3: (color online) Phonon dispersion relations of 
graphene (panel a), C- (panel b), B- (panel c) and W- (panel 
d) graphane. Acoustic and optical modes correspond, respec- 
tively, to blue (dark gray) and yellow (light gray) dispersions. 
Longitudinal and transverse acoustic branches are indicated 
as LA and TA, respectively. The acoustic branch with dis- 
placement patterns along the z-direction of Fig. [2] is marked 
as ZA. 
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while neither B- nor W-graphane are so. 

Finally, according to the present first principles total 
energy calculations we identified C-graphane as the most 
energetically favorable conformer. W- and B-graphane 
have higher ground-state energy of 0.05 and 0.10 eV (per 
C-H unit), respectively. These small differences in energy 
demonstrate that all three conformers are thermodynam- 
ically accessible, as indeed experimentally guessed^ 



IV. LINEAR ELASTICITY 

While C-graphane has trigonal symmetry (and, there- 
fore, is elastically isotropic as hexagonal graphene), the 
remaining B- and W-conformers show an orthorhombic 
symmetry, which causes an anisotropic linear elastic be- 
havior. Accordingly, the elastic energy density (per unit 
of area) accumulated upon strain can be expressed as 15 
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for the anisotropic ones. In Eqs.(TT]) and © we have ex- 
plicitly made use of the elastic linear constants Cn, C22, 
C12 and C44. Furthermore, the infinitesimal strain ten- 
sor e = |(Vu + VS T ) is represented by a symmetric 
matrix with elements e xx — e yy = and e xy = 

\ (jjf + ^7-) , where the functions u x (x, y) and u y (x, y) 

correspond to the planar displacement u = (u x , u y ). It is 
important to remark that Utrigo can be obtained from 
the U or tho by simply imposing the isotropy condition 
C11 = C22 and the Cauchy relation 2C44 = G\\ — C12, 
holding for both the hexagonal and trigonal symme- 
try. We will take profit of this by focussing just on the 
elastic behavior of a system described by Eq.©; when 
needed, the general results so obtained will be applied 
to the isotropic structures by fully exploiting the above 
conditions. The constitutive in-plane stress-strain equa- 
tions are straightforwardly derived from Eq.@ through 
T = dU /de, where T is the Cauchy stress tensor—. They 
are: T xx — C\\e xx + Ci2e yy , T yy = C 2 2^yy + C\ 2 (- X x and 

T X y = 2C44C xy . 

We now suppose to apply an axial tension a to any 
two dimensional hydrocarbon shown in Fig. [5] along the 
arbitrary direction n = cos9e x + siri9e y , where e x and e y 
are, respectively, the unit vectors along the zigzag and the 
armchair directions of the underlying honeycomb lattice. 
In this notation, therefore, 9 is the angle between ft and 
the zigzag direction. Under this assumption we get T = 
anigiri, where the in-plane stress components are defined, 



respectively, as T x 



a cos 2 9, TV 



a cos 9 sin 9, and 



T yy = a sin 2 9. By inverting the constitutive equation we 
find the corresponding strain tensor e. In particular, we 



easily get its longitudinal component e; 
the direction n 



n ■ e n along 
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as well as its transverse component et = t • e t along the 
direction t = — sin^e^ + cos9e y (with t ■ n = 0) 
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where A = C11C22 — C( 2 , c = cos6>, and s = sin#. By 
means of Eqs.(3) and (4) we obtain, respectively, the 
n-dependent Young modulus Eft = a/ei (i.e. the ratio 
between the applied traction and the longitudinal exten- 
sion) as 
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and the n-dependent Poisson ratio = — e t /e/ (i.e. the 
ratio between the lateral contraction and the longitudinal 
extension) as 



12(C 4 + S 4 ) 



Cn+C 22 -^)cV-C 
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Eqs.© and ([5]) are central to our investigation. 

First of all, we remark that they allow for the full linear 
clastic characterization of both the anisotropic graphane 
conformers and the trigonal one (as well as graphene), 
provided that in the latter case the isotropy and Cauchy 
conditions are duly exploited. In this case we in fact 
obtain the Young modulus E — (C^ — C 2 2 )/Cn and the 
Poisson ratio v = C12/C11, which are independent of the 
angle 9, confirming the planar isotropy. 

More importantly, however, Eqs.© and © imply that 
Efi and can be directly obtained by the linear elastic 
constants Cij , in turn computed through energy- vs-strain 
curves corresponding to suitable homogeneous in-plane 
deformations. This implies that there is no actual need 
to mimic by a computer simulation a traction experiment 
along the arbitrary direction identified by n or 9, indeed 
a technically complicated issue to accomplish. Rather, 
for the isotropic case (graphene and C-graphane) only 
two in-plane deformations should be applied in order to 
obtain all the relevant elastic constants, namely: (i) an 
axial deformation along the zigzag direction; and (ii) an 
hydrostatic planar deformation. For the anisotropic case 
(B- and W-graphane) two more in-plane deformations 
must be applied: (iii) an axial deformation along the 
armchair direction; and (iv) a shear deformation. The 
strain tensors corresponding to deformations (i)-(iv) de- 
pend by a unique scalar strain parameter £ as shown in 
Table HU For all imposed deformations the elastic energy 
of strained structures can be written in terms of Q as 



(7) 
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Table II: Deformations and corresponding strain tensors applied to compute the elastic constants Cj of graphane. The relation 
between such constants and the fitting term C/' 2 ' of Eq.(j7]l is reported as well. Deformations (i)-(ii) are applied to the C- 
conformer, while the full set (i)-(iv) of deformations is applied to the B- and W-conformers. £ is the scalar strain parameter. 

strain tensor [7 (2) U {2) 

isotropic structures anisotropic structures 

(oo) Cn Cl1 

2(C*n + C 12 ) C u + C22 + 2C12 

(°C) C ™ 

(° c $) 4C 44 



(i) zigzag axial deformation 

(ii) hydrostatic planar deformation 

(iii) armchair axial deformation 

(iv) shear deformation 



where Uq is the energy of the unstrained configuration. 
Since the expansion coefficient is related to the elas- 
tic moduli as summarized in Table [TTl a straightforward 
fit of Eq. has provided the full set of linear moduli for 
all structures. 

The synopsis of the calculated elastic constants for all 
graphane conformers, as well as graphene, is reported in 
Table IIII1 from which three qualitative information can 
be extracted. First, we observe that the difference be- 
tween Cn and C22 is much smaller for the B-conformcr 
than for W-graphane; therefore, this latter is by far the 
most elastically anisotropic conformer. Then, the value 
of C44, measuring the resistance to a shear deformation, 
decreases monotonically from graphene to W-graphane. 
Finally, we remark that the value of C\i (or, similarly, 
of the Poisson ratio) is much smaller in any graphane 
structure than in pristine graphene. The change in hy- 
bridization has therefore largely reduced the property of 
lateral contraction upon extension. Interestingly enough, 
the B-conformer is characterized by a negative C12 value, 
something unexpected and worthy of further investiga- 
tion, as reported below. 

Through Eqs.([5]) and © and by using the elastic con- 
stants reported in Table IIII1 we can quantify the in- 
dependence of E and v for the anisotropic structures 
by using polar coordinates, as illustrated in Fig. [4] and 
Fig. [5J respectively. In such a representation, a fully 
isotropic elastic behavior is represented by a perfectly 
circular shape of the £W and vn plots. This is indeed 
the expected, of graphene and C-graphane. On 

the other hand, Fig. @] confirms that W-graphane is much 
more anisotropic than the B-conformer. Furthermore, as 
anticipated, Fig. [5] provides evidence that the Poisson ra- 
tio in any graphane conformer is much smaller than in 



Table III: Graphene and graphane independent elastic con- 
stants (units of Nm -1 ). For graphene and C-graphane Cn = 
C22 and 2C44 = Cn — G 1 2 • 

graphene C-graphane B-graphane W-graphane 

Cn 354 248 258 280 

C22 225 121 

Ci 2 60 20 -1.7 14 

C44 93 81 



pristine graphene, since the corresponding va polar plots 
are contained within the graphene circle. 

An intriguing unconventional behavior is observed in 
Fig. \5\ for B-graphane, namely: for extensions along to 
the zigzag and armchair directions, the corresponding 
Poisson value is vanishingly small. This feature appears 
as a flower petal structure of the plot for such a sys- 
tem. By considering Fig. [SJ where a zoom of the pre- 
vious plot nearby the origin has been reported, we can 
actually learn more information. It is evident that four 
small lobes appear along the zigzag and armchair direc- 
tions (i.e. along the principal axis of the orthorhombic 
symmetry), corresponding to a Poisson ratio varying in 
the range —0.0075 < v < —0.0065. The limiting values 
are computed for extensions along the zigzag and arm- 
chair directions, respectively. It is truly remarkable that 
v could be negative in B-graphanc. While a negative 
Poisson ratio value is allowed by thermo-elasticity, this 
peculiar situation is only observed in special systems (i.e. 
foams, molecular networks or tailored engineering struc- 
tures) or just rarely in ordinary bulk materials (i.e. Si02, 
cubic metals, or polymer networks)^. 



o 




Figure 4: (color online) Polar diagram for the Young mod- 
ulus E of graphene and graphane conformers. The angle 8 
identifies the extension direction with respect to the zigzag 
one. Isotropic (anisotropic) behavior is associated to a circu- 
lar (non circular) shape of the En plot. 



6 



armchair 




zigzag 



graphene 
B-graphane 
C-graphane — 
W-graphane - ■ 



Table IV: Strain fields applied to compute the linear (Cij) 
and nonlinear (Cijk) elastic constants of the C-graphane. The 
relation between such constants and the fitting terms f/' 2 ' and 
[7 (3) of Eq.fTT]) is reported as well. 



Strain U (2) 




tensor 





(SS) 

(g°c) 

(So C ) 

(Vo) 



Cn Cm 

2 (Cii + Ci 2 ) 2Cm + 6C112 

2 (Cn — C12) 8C444 

3Cn — 2C12 Cm + I2C144 + 6C114 + 8C444 

3Cn — 2C12 -Cm — I2C144 + 6C114 + 8C444 

4 (Cn — C12) 12Cn4 — I2C124 + 8C444 



Figure 5: (color online) Polar diagram for the Poisson ra- 
tio v of graphene and graphane conformers. The angle 9 
identifies the extension direction with respect to the zigzag 
one. Isotropic (anisotropic) behavior is associated to a circu- 
lar (non circular) shape of the vn plot. The special case of 
B-graphane is enlighten by shading (see text). 



V. NONLINEAR ELASTICITY 



In this Section we generalize the previous analysis in 
order to draw a comparison between the nonlinear elas- 
tic behavior of graphene^ and the three conformers of 
graphane. The nonlinear strain energy function Uhex for 




zigzag 



Figure 6: (color online) The same as Fig. [5] zoomed in the 
region nearby the origin. Positive and negative Poisson ratio 
values are differently shaded as indicated. 



an hexagonal two dimensional lattice is^ 



V , 3 
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-Cii2el x e y 
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+ gClll4a 
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(8) 



where all the nonlinear features are described by the three 
independent moduli Cm, C222 and Cm- 

Similarly, the strain energy function Utrigo for C- 
graphane depending on the linear (Cn and C12) and 
nonlinear (Cm, C112, C144, Cm, C124 and C444) elastic 



Table V: Strain fields applied to compute the linear (Cij) and 
nonlinear (Cijk) elastic constants of the B- and W-graphane. 
The relation between such constants and the fitting terms 
U (2) and U {3) of Ea.ifTT]) is reported as well. 
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tensor 
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constants is found to be 

Utrigo = -Cii (el x + (? yy + 2el y ) + C12 (exxtyy ~ Z X y) 

+ g^lll ( e xx + e yy) + 2^ 112 ( e xx e yy + e xx^ yy ) 

+ 2Ci44 {(-xx^cy + e yy e xy) + ^114 ( e xx € xy + <? yy txy) 

+ 2Ci24 + -C^A^xy (9) 

For such a trigonal symmetry we have Cm = C 2 2 2) 
C\vi = C\22 and C144 = C244. Nevertheless, it is im- 
portant to underline that the overall nonlinear elastic 
response is truly anisotropic since not all the relevant 
isotropic conditions are fulfilled. 

Finally, the strain energy function U or tho for the B- and 
W-graphane, expressed through the linear (Cn, C22, G\i 
and C44) and nonlinear (Cm, C222, Cm, C122, C144 and 
C244) elastic constants, is given by 



Uortho = -jp\\t xx + -C 2 2tyy + 2C 4 4£ xy + Ci2€xx^yy 

+ gCme^ + -C 2 22£yy + 2 Cll2€ xx e yy 

+ -jC\22txx<?y V + 2C 14 ,4e xx e xy + 2C 2 44e vv e|j / (10) 

Eqs.©, © and (TTU)) can be obtained by using the stan- 
dard tables of the tensor symmetries, found in many crys- 
tallography textbooks (see for instance Ref. [TBh . 

As above described, in any symmetry the strain energy 
function depends on the third-order elastic constants (as 
well as the linear ones). Once again, they can be com- 
puted through energy-vs-strain curves corresponding to 
suitable homogeneous in-plane deformations. For each 
deformation the elastic energy of strained graphene or 
graphane can be written in terms of just the single de- 
formation parameter £ 

u(C) = u + \u^c 2 + |c/ (3) C 3 + o(c 4 ) (11) 

Since the expansion coefficients and are related 
to elastic constants, as summarized in Table IIVI for the 



Table VI: Graphene and graphane independent nonlinear 
elastic constants (units of Nm -1 ). 
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C- graphane 


B- graphane 


W-graphane 


Cm 


-1910 ± 11 


-1385±18 


-1609±31 


-1756±33 
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-1764 ± 3 




-1827±7 


-487±85 
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Figure 7: (color online) Polar representation of the nonlinear 
elastic moduli Dft of the three graphane conformers. In the 
B- and W-graphane cases, Da = D are everywhere negative 
(softening hyperelesticity), while in the C-graphene one the 
Da alternates negative and positive values (hardening hyper- 
elesticity). 



C-graphane and in Table IVl for the B- and W-graphane, 
a straightforward fit of Eq. (JTTJ) has provided the full set 
of third-order elastic constants. 

The results have been reported in Table WT\ where only 
the values of the independent elastic constants appear- 
ing in Eqs. (JU), (|9]) and (fT0|) are reported. We note 
that graphene and B-graphane are characterized by an 
inverted anisotropy: while Cm < C222 for graphene, 
we found C222 > Cm for B-graphane. On the con- 
trary, W-graphane has the same anisotropy of graphene 
(Cm < C222), but a larger |Cm— C222I difference. So, it 
is interesting to observe that the different distribution of 
hydrogen atoms can induce strong qualitative variations 
for the nonlinear elastic behavior of these structures. 

We finally observe that necessarily C444 = for B- 
and W- graphane because of the orthorhombic symme- 
try. On the other hand, this nonlinear shear modulus 
could assume any value for the trigonal lattice. Interest- 
ing enough, we have verified that C444 = also for C- 
graphane. This is due to the additional (with respect to 
the trigonal symmetry) mirror symmetry of C-graphanc. 

Similarly to the case of graphene,^2I a nonlinear stress- 
strain relation can be derived for the three graphane con- 
formers: 

<rn = Ene H + D H e% (12) 

where Efi and Dfi are, respectively, the Young modu- 
lus and an effective nonlinear (third-order) elastic mod- 
ulus, along the arbitrary direction n, as defined in 
Sec. IIVI The nonlinear elastic modulus Z)^ ri9 °^ for 
the C-graphane (as well as for any trigonal 2D lattice) is 



8 



given by 



jj(trigo) 



(C 1U -3Cii 2 ) 
+ (l-^)(l + ^ 2 )C m 

+6cs (1 + 1/) (1 + ^ 2 ) Cm - 12cs (1 + v) vC 12 ± 
+3c 2 s 2 (1 - v) (1 + i/ 2 ) (-Cm + 4Ci44 + C112) 
+4c 3 s 3 (1 + v) (1 + ^ 2 ) (— 3Cix4 + 2C444 + 3Ci 24 ) 
+8c 3 s 3 (l + ^)^ (-6C114 + 5C 4 44 + 6C124) 1 (13) 



while the corresponding modulus ]j^° rtho ^ g_ anc j \y_ 



graphene is 

j-y(ortho) 



1 



2A 3 El 



Cm (C22C 2 - C 12 s 2 



+C222 (Cns 2 - C 12 c 2 f 
+3C 112 (C n s 2 - C 12 c 2 ) (C22C 2 - C 12 s 2 ) 2 
+3C 122 (C22C 2 - C 12 c 2 ) (C llS 2 - C 12 c 2 ) 2 
-3C 166 c 2 s 2 (C22C 2 - C 12 s 2 ) (A/C44) 2 
-3C 266 cV (C llS 2 - C 12 c 2 ) (A/C 44 ) 2 ](14) 

Since Cijk < 0, as shown in Tab lVIl jj^° rtho ^ are negative 
for any direction ( see Fig. [TJ, so both B- and W-graphane 
show an hyperelastic softening behavior. The trigonal C- 
graphane behaves in a very different way instead. Since 
the Cm and C144 are positive, the C-graphane can show 



an hyperelastic hardening behavior in the angular sectors 
5/12tt + fcTr < 9 < 1/12 + fcTr and 8/12tt + fcvr < 9 < 
10/12 + kn (fceZ). 



VI. CONCLUSIONS 

In conclusion, present first principles calculations pre- 
dict that the class of auxetic materials is larger than re- 
ported so far, including as well two dimensional hydro- 
carbons like B-graphane. More precisely, since a negative 
Poisson ratio is observed for extensions along the zigzag 
and armchair principal directions, B-graphane is better 
referred to as an axially auxetic atomic sheet. Moreover, 
we calculated that the other two conformers, namely the 
C- and W-graphane, exhibit a vanishingly small value 
of the Poisson ratio. The nonlinear elastic behavior of 
graphane shows peculiar features as well. In particular, 
we have found that the C-graphane admits both soft- 
ening and hardening hyperelasticity, depending on the 
direction of the applied strain. These features makes 
graphane a very intriguing material with potentially large 
technological impact in nanomechanics. 
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